Can sterile neutrinos be ruled out as warm dark matter candidates? 
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We present constraints on the mass of warm dark matter (WDM) particles from a combined 
analysis of the matter power spectrum inferred from the Sloan Digital Sky Survey Lyman-a flux 
power spectrum at 2.2 < z < 4.2, cosmic microwave background data, and the galaxy power 
spectrum. We obtain a lower limit of ttiwdm ~ 10 keV (2a) if the WDM consists of sterile 
neutrinos and rawDM ~ 2 keV (2a) for early decoupled thermal relics. If we combine this bound 
with the constraint derived from x-ray flux observations in the Coma cluster, we find that the allowed 
sterile neutrino mass is ~ 10 keV (in the standard production scenario). Adding constraints based 
on x-ray fluxes from the Andromeda galaxy, we find that dark matter particles cannot be sterile 
neutrinos, unless they are produced by a nonstandard mechanism (resonant oscillations, coupling 
with the inflaton) or get diluted by some large entropy release. 

PACS numbers: 98.80.Cq 



Introduction. Warm dark matter (WDM) has been 
advocated in order to solve some apparent problems of 
standard cold dark matter (CDM) scenarios at small 
scales (see Q and references therein) , namely: the excess 
of galactic satellites, the cuspy and high density of galac- 
tic cores, the large number of galaxies filling voids. More- 
over, recent observational results suggest that the shape 
of the Milky Way halo is spherical |2( and cannot easily 
be reproduced in CDM models. All these problems would 
be alleviated if the dark matter (DM) is made of warm 
particles, whose effect would be to suppress structures 
below the Mpc scale. Detailed studies of the dynamics 
of the Fornax dwarf spheroidal galaxy suggest shallower 
cores than predicted by numerical simulations of CDM 
models and put an upper limit on the mass of a putative 
WDM particle [3| . One of the most promising WDM can- 
didate is a sterile (right-handed) neutrino with a mass in 
the keV range, which could explain the pulsar velocity 
kick 0, help in reionizing the universe at high redshift 
[j|, and emerge from manyparticle physics models with 
grand unification (e.g. 0,0)- Because of a small, non- 
zero mixing angle between active and sterile flavor states, 
X-ray flux observations can constrain the abundance and 
decay rate of such DM particles. The Lyman-a absorp- 
tion caused by neutral hydrogen in the spectra of distant 
quasars is a powerful tool for constraining the mass of 
a WDM particle since it probes the matter power spec- 
trum over a large range of redshifts down to small scales. 
In a previous work, |s| used the LUQAS sample of high 
resolution quasar absorption spectra to set a lower limit 
of 2 keV for the sterile neutrino mass. More recently, 
exploiting the small statistical errors and the large red- 
shift range of the SDSS Lyman-a forest data, Seljak et 
al. Isf found a lower limit of 14 keV. If the latter result is 



correct, a large fraction of the sterile neutrino parameter 
space can be ruled out (assuming that all the DM is made 
of sterile neutrinos); together with constraints from X-ray 
fluxes, this discards the possibility that DM consists of 
sterile neutrinos produced by non-resonant active-sterile 
neutrino oscillations (still, they could be produced by 
resonant oscillations caused by a large leptonic asymme- 
try in the early Universe , or considerably diluted 
by some large entropy release 0, 0, El, or generated 
in a radically different manner, e.g. from their coupling 
with the inflaton 01 )■ More recently, some joint anal- 
yses of the SDSS flux power spectrum and the WMAP 
year three data ^| have been presented in [Til . Hr| for 
standard ACDM models. The authors of ^ij found some 
moderate disagreement between the inferred power spec- 
trum amplitudes. Instead, from an independent analysis 
of the SDSS data [13, the authors of 0] find good agree- 
ment in their joint analysis. Here, we extend the analysis 
of ^3 to constrain the mass of WDM particles. 

Data sets and Method. We use here the SDSS 
Lyman-a forest data of McDonald et al. [13], which 
consist of 3035 quasar spectra with low resolution (R ~ 
2000) and low signal-to-noise spanning a wide range of 
redshifts (z — 2.2 — 4.2). The data set differs substan- 
tially from the LUQAS and C02 samples used || , which 
contain mainly high resolution, high signal-to-noise spec- 
tra at z ~ 2.5. More precisely, we use the 132 flux power 
spectrum measurements Pp(k,z) that span 11 redshift 
bins and 12 k— wavenumbers in the range 0.00141 < k 
(s/km)< 0.01778 (roughly corresponding to scales of 5- 
50 comoving Mpc). It is not straightforward to model 
the flux power spectrum of the Lyman-cv forest for given 
cosmological parameters, and accurate numerical simula- 
tions are required. McDonald et al. 0] modelled the flux 
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power spectrum u sing a large number of Hydro Particle 
Mesh simulations [18j, calibrated with a few small-box- 
size full hydrodynamical simulations. Here, instead, we 
model the flux power spectrum using a Taylor expansion 
around a best fitting model: this allows a reasonably ac- 
curate prediction of the flux power spectrum for a large 
range of parameters, based on a moderate number of full 
hydrodynamical simulations |lf)j . The method has been 
first introduced in Ref. [lij and we refer to this work 
for further details. The fiducial flux power spectrum has 
been extracted from simulations of 60 h~ 1 comoving Mpc 
and 2 x 400 3 gas and DM particles (gravitational soften- 
ing 2.5 h" 1 kpc) corrected for box size and resolution 
effects. We performed a number of additional hydrody- 
namical simulations with a box size of 20 hr 1 comoving 
Mpc and with 2 x 256 3 gas and DM particles (grav. soft. 
1 h^ 1 kpc) for a WDM model with a sterile neutrino 
of mass m s = 1,4,6.5 keV, to calculate the flux power 
spectrum with respect to changes of the WDM particle 
mass. We have checked the convergence of the flux power 
spectrum on the scales of interests using additional sim- 
ulations with 2 x 256 3 gas and DM particles and box 
sizes of 10 h^ 1 Mpc (grav. soft. 0.5 hr 1 kpc). We then 
used a modified version of the code CosmoMC [2lj to 
derive the parameter likelihoods from the combination of 
the Lyman-a data with Cosmic Microwave Background 
(CMB) and galaxy power spectrum data, from WMAP 
Q, ACBAR H3, CBI H3, VSA [H and 2dF 0. In 
total, we used a set of 29 parameters: 7 cosmological 
parameters; 1 parameter describing a free light-to-mass 
bias for the 2dF galaxy power spectrum; 6 parameters 
describing the thermal state of the Intergalactic Medium 
(paramctrization of the gas temperature-gas density re- 
lation T — To(z)(l + (5) 7 W -1 as a broken power law at 
z — 3 in the two astrophysical parameters Tq(z) and 
7(2)); 2 parameters describing the evolution of the ef- 
fective optical depth with redshift (slope and amplitude 
at z = 3); 1 parameter which accounts for the contri- 
bution of damped Lyman-a systems and 12 parameters 
modelling the resolution and the noise properties (see 
|2fij). We applied moderate priors to the thermal history 
to mimic the observed thermal evolution as in |2fij . but 
the final results in terms of sterile neutrino mass are not 
affected by this. 

Results. Wc assume the Universe to be flat, with 
no tensor or neutrino mass contributions. We further 
note that adding CMB and large scale structure data 
has very little effect on the results for m s , since the free- 
streaming effect of WDM particles is visible only on the 
scales probed by the Lyman-a flux power spectrum |j3(| . 

In Figure ^ we show the 2-dimensional marginalized 
likelihoods for the most important cosmological and as- 
trophysical parameters: as, n s , S7 m and the effective op- 
tical depth amplitude measured at z = 3, T^g, all plotted 
as a function of the parameter (1 keV)/m s . The con- 
straints on m s get stronger for the Lyman-a forest data 
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FIG. 1: 2-dimensional marginalized likelihoods (68% and 
95% confidence limits) for n 3 ,as, fi m and the effective optical 
depth at z = 3, using the SDSS data at z < 4.2 (left, green), 
z < 3.6 (middle, white) and z < 3.2 (right, blue). 



in the highest redshift bins. To demonstrate this we plot 
the likelihood contours for data in three different red- 
shift ranges: z < 3.2 (blue), z < 3.6 (white) and z < 4.2 
(green), which is the whole data set. The constraints 
improve by a factor almost 3 (2) for the whole data set 
compared to the z < 3.2 (z < 3.6) subsamples. At high 
redshifts, the mean flux level is lower and the flux power 
spectrum is closer to the linear prediction making the 
SDSS data points very sensitive to the free-streaming ef- 
fect of WDM 0. We find no strong degeneracies between 
m s and the other parameters, showing that the signature 
of a WDM particle in the Lyman-a flux power is very 
distinct, and that other considered cosmological and as- 
trophysical parameters cannot mimic its effect. 

In Figure [5] wc show the 1-dimensional marginalized 
likelihoods for (1 keV)/m s for several redshift ranges. 
The 2a lower limits for the sterile neutrino mass are: 3.9, 
8.3, 8.1, 8.6, 10.3 keV for z < 3.2, 3.4, 3.6, 3.8, 4.2, respec- 
tively. The corresponding limits for an early decoupled 
thermal relic are: 0.9, 1.7, 1.6, 1.7, 1.9 keV (see @| for 
the correspondence between the two cases). Also shown 
(dotted black line) is the constraint obtained in || using 
the LUQAS and C02 samples The SDSS data 

improve the constraint from the high resolution data at 
z ~ 2.5 by a factor 5. This is mainly due to the extension 
to higher redshift where the flux power spectrum is most 
sensitive to the effect of WDM. The smaller statistical 
errors of the flux power spectrum and the coverage of a 
substantial range in redshift help to break some of the 
degeneracies between astrophysical and cosmological pa- 
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FIG. 2: 1-dimensional marginalized likelihoods for the pa- 
rameter (1 keV)/m s for the SDSS Lyman-a data for the red- 
shift ranges z < 3.2, 3.6, 3.8, 4.2 and the VHS [gjj data. 



rameters and also contribute to the improvement. Our 
independent analysis confirms the limits found in 9] for 
the SDSS Lyman-a data and a small sample of high res- 
olution data that also extends to high rcdshift. Note, 
however, that our lower limit for essentially the same 
data set is ~30% smaller (indeed, when using only SDSS 
Lyman-a data, Ref. obtains m s > 12 keV (2a), which 
includes a 10% correction caused by the non-thermal mo- 
mentum distribution of sterile neutrinos |28j : so, for the 
assumption made here, they would get m s > 13 keV). 

Discussion. In Figure [3] we summarize a num- 
ber of current constraints for sterile neutrinos in the 
(m s , sin 2 2(9) plane, where 9 is the vacuum 2x2 mix- 
ing angle between active and sterile neutrinos [3(|. We 
show the limits obtained from different types of X-ray 
observations: X-ray diffuse background (XRB, o rang e 
curve, [HI); flux from the Coma cluster (blue curve, [12); 
and finally, flux from the Andromeda galaxy (M31) halo 
(95% C.L., green dashed curve, H^). In addition, we 
plot the Lyman-a constraints obtained in this work (red 
dashed) and in [jj (black dotted). The region which 
can explain observed pulsar kicks M is shown as the 
hatched area. Finally, according to y|, sterile neutrinos 
produced from non-resonant oscillations (i.e., in absence 
of significant leptonic asymmetry, L = 0) with a den- 
sity Odm = 0.23 ± 0.04 should lie between the two black 
solid curves (the computation in is based on simpli- 
fying assumptions concerning the QCD phase transition; 
the effect of hadronic corrections is currently under in- 
vestigation [34) and could shift the allowed region in the 
(m s , sin 2 29) plane). If all these constraints are correct, 
then there is no room for sterile neutrinos as DM candi- 
dates in the standard case. Models in which the decay 
of massive particles release some entropy and dilutes the 
dark matter by a factor S can alleviate the tension be- 
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FIG. 3: This plot summarizes some of the parameter space 
constraints (at the 95% C.L.) for the sterile neutrino models, 
assuming that they constitute the dark matter. Limits are 
explained in the text. 



tween the Lyman-a and X-ray bounds ^lj; but a very 
large S is needed 0, 0] . As mentioned in the introduc- 
tion, the sterile neutrino remains a viable WDM candi- 
date for alternative production mechanisms (e.g. reso- 
nant oscillations with L ^ 0, or coupling with the infla- 
ton). Recently, Ref. 01 questioned the results based on 
the LMC and Milky Way because of uncertainties in mod- 
elling the dark matter distribution; and also those based 
on detecting emission lines in cluster spectra |32j , which 
used a fixed phenomenological model for X-ray emission 
(no t shown in the figure but 30% more constraining than 
|3l|). If these observational constraints are inaccurate, 
then a sterile neutrino mass in the range 955 m s (keV) 
11.5 and sin 2 29 ~ 2 x 10 -9 would be marginally consis- 
tent with the XRB bound and the Lyman-a forest data, 
but it is strongly excluded by the robust limit obtained 
by Ref. [3^ (which is very conservative, since the bound 
quoted as 2a by the authors requires a signal a few times 
larger than the background). The corresponding emis- 
sion line for such a decaying sterile neutrino would be at 
E ~ 5.5 keV (close to, or possibly contaminated by, the 
recently discovered Chromium line (35| ). If instead all 
X-ray constraints are correct, but the two recent Lyman- 
a forest constraints are not accurate, then a mass of 
m s ~ 2 keV is still possible and compatible with the 
robust and conservative lower limit from ||. It would 
also satisfy the requirement from the dynamical analy- 
sis of the Fornax dwarf galaxies Q. However, the latter 
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possibility appears unlikely. Even if the highest rcdshift 
bins of the SDSS Lyman-a forest data were affected by 
not yet considered systematic errors the analysis of the 
data with z < 3.2 still gives a lower limit of about ~ 3.5 
keV (see HJ). Appealing to an insufficient resolution 
of the hydrodynamical simulations would also not help, 
since an increase in resolution could only increase the flux 
power spectrum at small scales and raise the lower limits. 
We have furthermore checked explicitly that this is not 
the case and that other possible effects on the flux power 
have a different signature than that of WDM. A poten- 
tially big improvement on the quality of the constraints 
from Lyman-a forest data could be achieved by an anal- 
ysis of a large set of high-redshift, high-resolution data to 
extend the measurement of the flux power spectrum at 
high redshift to smaller scales. This would, however, also 
require accurate modelling of the thermal history and the 
contribution of associated metal absorption to the small 
scale flux power spectrum. 
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